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We apply the recently proposed susy Hubbard operators [P. Coleman, C. Pepin and J. Hopkinson, 
Phys. Rev. B 63, 140411(R) (2001)] to an atomic model. In the limiting case of free spins, we 
derive exact results for the entropy which are compared with a mean field + gaussian corrections 
description. We show how these results can be extended to the case of charge fluctuations and 
calculate exact results for the partition function, free energy and heat capacity of an atomic model 
' for some simple examples. Wavefunctions of possible states are listed. We compare the accuracy of 

large N expansions of the susy spin operators [P. Coleman, C. Pepin and A. M. Tsvelik, Phys. Rev. 
B 62, 3852 (2000); Nucl. Phys. B 586, 641 (2000)] with those obtained using 'Schwinger bosons' 
and 'Abrikosov pseudo-fermions'. For the atomic model, we compare results of slave boson, slave 
£f>^ , fermion and susy Hubbard operator approximations in the physically interesting but uncontrolled 

limiting case of N— >2. For a mixed representation of spins we estimate the accuracy of large N 
expansions of the atomic model. In the single box limit, we find that the lowest energy susy saddle- 
. point reduces to simply either slave bosons or slave fermions, while for higher boxes this is not the 

case. The highest energy saddle-point solution has the interesting feature that it admits a small 
region of a mixed representation, which bears a superficial resemblance to that seen experimentally 
close to an antiferromagnetic quantum critical point. 
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' ■ 1. INTRODUCTION 

One of the interesting challenges arising in systems exhibiting strong electronic correlations is the ability to ade- 
quately develop a microscopic description of metallic physics in. regions where Fermi liquid theory breaks down. The 
observations of & scalingcl, linear resistivity over three decadesa in temperature coinciding with a Tln[-pr] divergence 
| in the specific heat capacity!!! in YbRIi2Si2 is reminiscent of properties shown by CeCug-ajAua; which additionally 
. exhibits ^ scalingtil. In these heavy fermion systems close to an antiferromagnetic quantum critical point only one 
energy or length scale seems to be important-which generally only occurs in theories which are below their upper 
critical dimensions Evidence for a qualitatively similar phase diagram for the high-T c cuprates has been advanced 
by Tallon and LoramQ who have assembled a large amount of experimental data on T* and interpret it in terms of 
Ji-j an underlying quantum critical point, presumably responsible for the linear T dependence, seen in the normal state 
■ of these systems that has given rise to. the marginal Fermi liquid phenomenology of VarmaQ and others. Recent /LiSR 
measurements by Panagopoulos et ajQ and an as yet unsubstantiated report of high temperature superconductivity 
"j~>; . in a FET of CaCu02 by Schon et al3 have supported this interpretation. Again, the culprit for this unusual metallic 
behaviour is believed to be an underlying quantum critical point. If the entire Fermi surface fundamentally breaks 
down in this region, one might imagine that the competition between spins trying to magnetically order and those 
trying to form heavy quasiparticles could lead to a new kind of quasi-particle excitations That is, if the Fermi tem- 
perature approaches as we approach a quantum critical point, we may see the formation of novel kinds of states; 
perhaps states of higher symmetry classes could become important even at these low temperatures. The key question 
that must be addressed is: how do spin and charge interact at the brink of magnetism? 

Theoretical approaches to this problem are hindered by the difficulty of capturing the profound transformation in 
spin correlations that develops at the boundary between antiferromagrijetism and paramagnetism. Usually we model 
this behavior by representing the spin as a boson in a magnetic phaseja or as a fermion in a paramagnetic phase,c2l 
but by making this choice, the character of spin and charge excitations which appear in an approximate field theory is 
restricted and lacks the flexibility to describe the co-existence of strong magnetic correlations within a paramagnetic 
phase. 

These considerations have motivated the development of new methods to describe the spin and charge excitations 
of a strongly correlated material which avoid making the choice between a bosonic or fermionic spin.liiJiiil Recently, we 
introduced a new representation of the Hubbard operatorsEj which attempts to treat magnetism and paramagnetism 
on an equal footing within an approximate field theory allowing a description of physics in terms of operators which 
obey canonical commutation relations, yet avoid forcing us to specify the nature of a system's ground state. In this 
approach, we introduce vector fields 

F a = (h,-,fN,4>) B a = (bi,...,b N ,x) (1) 
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where f a and b a are-Abrikosov pseudo-fermionsQ and Schwinger bosonsi respectively—while 4> an d X are slave bosonsEzI 
and slave fermionsEJ. For \a> £ \a >, |0 >, the Hubbard operators can be written!!! 

\a >< b\ = X ab = B\B b + F\F h (2) 

To guarantee that the Hubbard operator representation is irreducible, we need to set the values of the linear and 
quadratic Casimirs of the group (here g — Diag[l ... 1, —1]), 



C* (1) = Tr[X] = Q, C {2) = Ti[XgX] = Q(N — 1—Y) 



(3) 



The first constraint corresponds to summing the total number of spin and hole states at a site (henceforth labeled by 
Q), while the latter describes the degree of antisymmetry of a given representation (which is fixed by Y = h-w (see 
Fig.l)). In terms of canonical operators, we can express Q and Y as: 



= n b + rt0 + rif + n x , 



(4) 



Y 



+ 71/ - (n b + n x ) + — [0,0 f ], 



(5) 



where = b\f a — x' is an operator interconverting fermions and bosons for the corner state. While the physically 
interesting case (Fig.l (a)) corresponds to Q=l, Y=0 and N=2, one can hope to learn new qualitative features from 
large N expansions about different symmetry classes as has been done previously for Q = Qo, Y = 1-Qo ((Fig.l(c))slave 
fermions-magnetism) ; and Q = Qo, Y = Qo - 1 ((Fig.l(d))slave bosons-paramagnetism) . One now can also treat 
cases like (Fig. 1(b) Y=0 -paramagnetism + magnetism?). 
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FIG. 1. (a) Fundamental representation (Q, Y) — (1,0), (b) L-shaped Young tableau corresponding to the spin representa- 
tion generated by supersymmetric Hubbard operators. The asymmetry Y = h — w and Q is the number of boxes, (c) Young 
tableau for fully symmetric representation corresponding to Slave fermion limit (d) Fully antisymmetric, slave boson limit. 



In this paper we first consider the limit of free spins with strong correlations in which charge fluctuations are totally 
suppressed. In this limit, Q becomes the total number of spins (n b + nf) and the energy of the empty state is taken 
to infinity. We derive exact results for the entropy of a general Q,Y. Using a mean field approach plus gaussian 
corrections, we compute the entropy of the general free spin case and compare this with -^expansions of the exact 
results (Stirling's approximation). We then compare the supersymmetric spin formalipsKpEJ with earlier large N 
approximations. Representations of spins in terms of N-component 'Schwinger' bosonsEilcl and 'Abrikosov pseudo- 
fermions'tJ become exact in the large N limit, but in the physical limit N — * 2, provide only approximate descriptions 
of the physics. Calculation of the exact entropy of free spin states, allows us to estimate the error inherent in mean 
field descriptions, mean field plus gaussian corrections, and self-consistent mean field + correction approximations to 
each of the three methods. Surprisingly, self-consistent approximations seem to not always do as well as mean field 
approximations at low N. 

Re-introducing charge fluctuations, we next treat an atomic model. In such a model there is an excitation energy for 
a spin to leave its localized site, creating an empty or vacuum state. The Hamiltonian for such a system is particularly 
simple, being given by simply the energy difference between the unoccupied site and the occupied spinful state. In 
terms of Hubbard operators, this can be written: 
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H — EdX acr — EdQ — EdX 00 — —E^Xoo 



(6) 



where Ed is the energy of the localized level, and we have dropped a constant term in the last term for convenience. 
Now, we have charge fluctuations in addition to spin fluctuations allowed by the system. For small N we can write 
exact wavefunctions corresponding to possible states of the system and their associated energies. Then, we can write 
correspondingly exact free energies, entropies and heat capacities for such a system. These exact values allow us at 
low N to compare the efficacy roi the supersymmetric Hubbard operator formalism with that of previously existing 
slave bosonlHI and slave fermiorJIS approximations. We find that even in this uncontrolled regime, bosonic descriptions 
(slave fermions) provide the best approximation to spin states, while fermionic descriptions of the spin (slave bosons) 
provide the best approximation to charge states at N=2. The lowest energy saddle-point of the susy solution at N=2 
corresponds to this solution, while the highest energy saddle-point generates a phase diagram which bears a superficial 
resemblance to the phase diagram of the cuprates and includes a small region where the solution is neither fermionic 
nor bosonic but mixed. For higher N, we show how one can generate the exact partition function and estimate the 
accuracy in calculations of the entropy and heat capacity of this new large N approach which allows representations 
of the spin neither bosonic nor fermionic in character. 



2. FREE SPIN LIMIT 



First, let us consider the free spin. This corresponds to taking the limiting case when we send Ed — > — oo. 



A. Entropy 
1. Exact Result 



We are left with SU(N) spins as the only states of our system. To count the number of states-available, we make 
use of the fact that the dimension of an irreducible representation of SU(N) can be expressed as:Ej 

A „= n <im^> (7) 

1,3 

where gij is the hook length of the (i,j) box of the Young tableau (see Fig. 2 a)) defined as the number of neighbors 
to the right and bottom of the box (see Fig. 2 b)). For the specific case of L-shaped Young tableaux, this yields, 
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FIG. 2. (a) An L-shaped Young tableau with w boxes horizontally, and h boxes vertically, (b) The hook length of each box, 
defined as the number of neighbors to the right and bottom of the box itself.E3 



From this formula, one can deduce an exact expression for the entropy depending only on N (SU(N)), w (width of 
the Young tableau), and h (height of Young tableau). This equation is: 



5 = In 



N \ / N + w 



h J \ w J N(w + h-l) 



(9) 



Since we work in the large N limit it would be useful to know the leading dependence of this exact result. We can use 
Stirling's approximation to expand the factorials, i.e. 



ln(d) = xln(-) + - \u(2ttx) (10) 
e 2 



which, after some work, implies that 



SstiHing = N(S f [h] + S b [w]) - \ (In {2-kNw{\ +«)-!)) + H2nNh(l - h))) + H N{w 'f h _ 1) ) (H) 

where h= ±,h= ^.w^ § ,w = ^±1 and 

Sf[x] = —[xhix + (1 — x) ln(l — x)] 
S b [x] = -[xhix- (1 + x) ln(l + x - jf)} (12) 

In order to express h, and w as functions of N, we can remove the leading dependence of Q and Y on these variables 
and expand in it. We can choose the new variable (other than N) in such a way that its sign agrees with that of 
Y, and its maximal amplitude is 1 for the special case when Q = f . That is, we define the variable £ such that 
w = ^(1 — £) + and h — ^(1 +£) + "2 ' With this assignment, for the case we will consider (Q — y), — 
and Y = (^ — 1)£. The physics at N = 2 will be invariant of the value of £. It does have the problem that N — > 00 
and £ — > ±1 do not commute. 

Following the above-developed procedure, we express the exact entropy of the description as (Q = y): 



2, f N\ f N + w-l\ wh 



(13) 



N N~" \\ h J \ w J N(f) J 
and expansion of this according to Stirling's approximation gives us: 

ZSstirling A , t- A lu, , 1 \ , t-A ^ \\ 

~^T~ = - ( 2 + % - iV ))ln(( 4 + 2N )+ ^4 - 2iV )} 

- (§ - 4 - ^)) Mil - 2A7) - «i - ^) - (5 + 4 - S» ln(( I + 2^ 

+ f(— - -)) + (- + £( — - -)) ln((- - — ) + 6(— --))-— ln(7r7V) - — ln(l - — ) (14) 
y 2N \" y 2 y N 2" u 4 2N' 2N A" N y ' N v N' ' 

extracting terms of order from inside the logs, this becomes the leading order approximation: 
2S SU r^ lead = + f(1 _ | )} ln(1 + _ 1 (3 + + | )} ln(3 _ ^ 

- \A + + |)) Ml - + + + |)) ln(5 - - I In(vr^V) - | (15) 



2. Mean field + Gaussian Corrections 



(a) Susy. Since we cannot enforce the constraints Q = Qo and Y = Yq rigorously at the mean field level, (Y has a 
four-fermion term of order jj), we would like to see how inclusion of these terms consistently affects the saddle point 
approximation we make. For free spins, inclusion of the constraints as Lagrange multipliers leads to an action of the 
form 
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S = f dr(J2 Udr + A + Q.U + £ h(dr + A - Ob a + £ ^{hUffo - f<,b a bpf ) - AQ - Ob (16) 

The total free energy then is made up of five terms to order A. , 

Ftot = J V[f, J, b, b]e- s = F m/t + F^ g . c . = AF/ + 7VF b + F v - \Q - (Y + F 5Xf + F 5Xb (17) 
where the last three terms correspond to calculating contributions due to 

SA, f 

, and U SXb = ^» J , (18) 

where the two terms of order N would ordinarily define the location of the saddle point, and the three latter terms 
be evaluated at the location of the mean field saddle point. As this would correspond to setting -j < [6, 9<] >= at 
the saddle point, we choose to include F v in determining the location of the saddle point itself. This has the effect of 
coupling the gaussian fluctuations in A/ and Af, such that the sum F$\ f + Fg\ b is replaced by a single term Fs x s Xb - 
In the above, A/ = A + £ and A b = A — £ arise from Lagrange multipliers A and £ which have been introduced to 
enforce Q = h+w — 1 and Y = h — w respectively at the mean field level. r\ is a Hubbard-Stratonovich field introduced 
to decouple the O(-j^) term from Y. Details of the calculation are given in Appendix A. Grouping all our terms from 
the Gaussian approximation to the entropy (calculated in Appendix A), we have 

F = -ATlnfl + e- 0X f) + NT ln(l - e"^) + Tln( nfe _ + . n/ ) - \(h + w) - ((h - w) 

2h b hf 

+ - ln((27r) 2 (A^ b (l + h b )Nh f (l - hf) - n a (l - n a ){Nn b (l + h b ) + Nn f {l - n/)))) (19) 



Since there is no energy term here F = —TS to calculate the entropy we simply need to divide by the temperature 

OF I _ _dF_ _ dF d\ I 
dT I" — dT d\ dT I 



although we can also perform the constrained derivative S — — §fJ|n = — fy — §xwU- 



- 2 ,n,l + „-*,, - 21n(l - e -) + + - 2 ln(=|^) 



- -£ ln((2^) 2 (An b (l + h b )Nh f {l - hf) - n„(l - n a )(Nh b (l + h b ) + Nh f (l ~ h f )))) + -|(ln(2)) (20) 

where we have used the extra —j3\ from the constraint on Q to create the third to last term as shown in Appendix 
A. Here we have n a = e p 2 l +1 , n f = 0>-) +1 an< ^ = e 3A t--i wnere -V = ^ + C an( l A b = A — ( are the Lagrange 
multipliers as they affect fermions and bosons respectively. The last term we remove as the factor of two from the 
double counting of the corner box. Since 2h — q + y + and 2w = q — y + A, the extra A contribution allows us to 

neatly group prefactors in terms of h and w. That is: 

= -2(1 - h) ln(l - hf) - 2h\nh f + 2(1 + w) ln(l + n b ) - 2wlnh b - — \n( nb + Uf ) 
N N h b hf 

- — ln((2^) 2 (An b (l + h b )Nh f (l -h f )- n a (l - n a ){Nh b {\ + h b ) + Nh f {l - h f )))) (21) 
Taking d ^ !t = = d ^ b ft where the mean field contribution to the free energy shown in Eq. (18) is given by: 

F mft = NT((1 - h) ln(l - h f ) + hln(h f ) - (1 + w) ln(l + h b ) +w\n(h b )) + T\n{ U ^ +nf ) (22) 

n b rif 

yields the saddle-point equations: 

h f + i(l -n a )=h (23) 
hb + jjn a = w (24) 
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where n a = e2 p\ +1 = zrr^n — 1S the gaussian contribution to the — ^ < [6,6^] > term of the constraints which we 



n b 

need to include in the saddle-point (ie. at the mean field level). With this assignment, we see that at the mean field 
level we recover the results given by the earlier Abrikosov pscudofcrmion and Schwinger boson approaches. Analytic 
solution of these equations shows a recovery of the constraints for these two cases when the Young tableau is a column 
or a row respectively, while for Y = we have in general a mixed case (see Fig. 3b). 

(b) Abrikosov fermions. The exact entropy (normalized) of a column of height ^ an d width 1, is given by: 



25 2 f N \ 

iV=A7 ln U) (25) 
Using Stirling's approximation to expand the factorials, we obtain 

2 ^7" 9 ' S =21n2-lln(^) (26) 

which should be the best that a jj theory can do. With self-consistent corrections, one might hope to approach the 
exact solution. 

For Abrikosov pseudofermions the mean field + gaussian corrections free energy is given by: 

F = —NT ln(l + er fJX ) - XQ + \ ln(27rJVfi/(l - ft/)) 

= iVTln(l - ft/) - NqT\n{^-) + | ln(27rAft/(l - ft/)) (27) 

where the last term is the gaussian corrections to the mean field result. Taking the usual mean field saddle-point 
9 $™p = (mfc) we obtain ft/ = q, such that the entropy (at q = |), 

^=21n2-lln(f ) (28) 

reproduces the Stirling's approximation to the exact result. An attempt at a fully 'self-consistent' (fsc) description 
„_l i 

might set ft r = 2 f but such attempts do not seem to improve the accuracy. 

(c) Schwinger Bosons. For the simple case of a purely bosonic description of spins, the exact entropy (for a row of 
width y and height 1) is given by: 

2S exact 2/Ar + |-l\ (29) 



N N 

which using Stirling's approximation simplifies to: 

= (3 - i) ln(| - i) + In 2 - j, \n(nN) + (-2 + i) ln(l - i) (30) 
which to leading order in N (ie dropping inside log's) yields 

2Sg "7 g '' ead = 3 In(|) + ln(2) - 1 In(^) (31) 
The mean field + gaussian corrections free energy is given by 

F = NT ln(l - e~ px ) - XQ + — ln(27rAft fc (l + n b )) 

= —NT ln(l + h b ) - NqT\n{^-^-) + - ln(27rAft b (l + n h )) (32) 

n b 2 

For the usual mean field saddle-point, we obtain hb = q so the entropy is found to be: 

S m fc = N(l + q) ln(l + q) - Nq ln(q) - 1 ln(27rNq(l + q)) (33) 
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which reduces to 

2S mfc Q1 .3. 1 .3ttJV 



31n(-)+ln2--ln(^-) (34) 



for the special case q = h, agreeing with the leading approximation to Stirling's approximation to the exact entropy. 

Q+ — 

For a fully 'self-consistent' saddle-point the consistency conditions would yield h — 2 f which is less close to the 

1 JV 

known exact result. 

B. Comparison of Accuracy 
1. Susy 



Susy spin: Entropy vs. 1/N for Y = 

h = N/2, w = N/2 non-physical divergence at N=2 



Occupation vs. N for h/N = w/N = 1/2 

Free Spins Analytic (N=2 admits also n/N=0,n b /N=0.5,n a =0) 
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FIG. 3. a) As N increases the supersymmetric spin mean field + gaussian corrections at Y = approaches the exact result, 
where the saddle-point has been chosen to maximize the number of terms kept at the mean field level. Choosing the incorrect 
saddle-point (ignoring the 1/N term in the constraint for Y) underestimates the entropy. About the correct saddle-point the 
divergence as one approaches N=2 is a non-physical term arising from the gaussian fluctuations in £ and A and might be expected 
to be cancelled by fluctuations in the corner box. Note that the correct saddle-point reproduces the Stirling's approximation 
when sufficiently far from the non-physical divergence at N=2. At N=2, Schwinger bosons and Abrikosov fermions solutions 
are found, but the entropy has an additional non-physical — ■x ln(47r0) term, b) In the limit N— >2, either Schwinger bosons or 
Abrikosov fermions are produced by the mean field conditions. For N>2 there is a unique analytic solution at Y=0 for h = 1/2 

At N=2, the susy approximation reverts to either Abrikosov fermions or Schwinger bosons (see Fig 3(b)) recovering 
the entropy of these approximations if one is willing to remove a divergent term by hand: — j^ln(2irN(0)). It is 

interesting to speculate that a correction of the form — -^log( ) ma y be necessary to remove this 

divergence in these fluctuations and might arise in higher order diagrams, but such a correction is beyond the scope 
of this paper. Nonetheless, the agreement for N>2 is quite reasonable (see Fig. 3(a)). 

2. Abrikosov fermions 



To compare these approximate solutions with the exact solution, by N — 10 a graph is indistinguishable by the eye 

T 



(see Fig. 4(a)). At JV" = 2, 2Se |° ct = In 2, so the percentage error is given by 



2.S' 



Stirling' s,mfc 



2S, 



1 — * 100 = 17% 



(35) 



which agrees with the fully self-consistent version for q = \ . So there seems to be no advantage here to fixing 9 g^ BC = 
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Schwinger bosons: Entropy vs. 1/N 
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FIG. 4. Beyond N~ 12, Abrikosov fermion and Schwinger boson solutions for the entropy of a free spin are indistinguishable 
to the eye from the exact for a free spin. Here we compare the exact entropy as a function of spin degeneracy N two large 
N approximations. For Abrikosov fermions, at the mean field saddle-point, mean field + gaussian corrections results agree 
with Stirling's approximation to the exact result and an attempt to introduce self-consistency. The slave boson mean field 
saddle-point out-performs both the full Stirling's approximation and an attempt at full self-consistency. 



3. Schwinger bosons 



The self-consistent approximation ( d g^ sc = 0) does dramatically worse (—40% error at N =2) than mean field 



solution = 0). At N = 2, 



2S mfc 25 c 



2 



2S e 



* 100 = 14% (36) 



whereas 

2 S SUrlin „, s 2S e 



2S e 



* 100 = 17% (37) 



while for large N, both approximations approach the correct result, above N = 12 being indistinguishable (see 
Fig.4(b)). It is interesting to note that the leading -k approximation to Stirling's formula does a better job than the 
formula itself, and the mfc approach reproduces this nice feature at N=2. 

3. ATOMIC MODEL 

Next, we consider a simple atomic model. In such a model there is an excitation energy for a spin to leave its 
localized site, creating an empty or vacuum site. The Hamiltonian for such a system is particularly simple, being 
given by simply the energy difference between the unoccupied site and the occupied spin-ful site. In terms of Hubbard 
operators, this can be written: 

H = —EdX o- (38) 



A. Entropies, Heat Capacities and Wavefunctions 
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1. Exact Results 



(a) N— 2. For a single SU(2) spin with no double occupancy, we know that the exact partition function is given by 

Z = 2 + e< 3Ed (39) 

which means that the free energy is given by 

F = -T\n{2 + e f)Ed ) (40) 

the entropy is 



and the heat capacity is given by 



_ 2(fiE d ) 2 e- (iEd 

v ~ (2 e -^ + l) 2 1 j 



By writing out the possible states of the system and their energies, we can count the total number of states available 
per site, as done for the Q = 1, Y = 0, N = 2 case above. This procedure allows us to compute exact results for 
thermodynamic quantities for larger representations. In Table I, we list the results for Q=2,3 calculated in Appendices 
B, C. 



Q,Y F S = -$ C v = -rffi 

{2,1} -Tln(l + 2e^+ e 2 ^) ln(l + 2e^ + e 2 ^) + 'gt'-Xn* 

{2,-1} -Tln(3 + 2e^) ^(3 + 2^ + 5^5 T^X^ 

{3,0} -Tln2-rin(l + 2e^+e 2 ^) In 2 + ln(l + 2e? E * + e 2 ^) + ^X^ 

{3,-2} -Tln(4 + 3e^) ln(4 + 3e^) + ^f»l + 7 

{3,2} -Tln(l+2e^ +e 2 ^) ln(l + 2e pE " + e 2 ^) + ^ggf^ ^ffX^T 

TABLE I. Free energy, entropy and heat capacity values of Q, Y restricted to the case when N = 2. The first non-trivial 
case with Y=0 has Q = 3. 
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(b) Exact Wave] 'unctions. For the simple case (Q=l, N=2) we can choose to represent the spin as either a fermion 
or a boson, as in reality a single spin is neither. This also applies to the empty state. It can be described by a 
grassman variable or a spinless boson. This dictates the possible wavefunctions for Q = 1, Y = 0, N = 2 to be: 

/||0>, /J|0>, &{|0>, &{|0>, f |O>, x + |0 > (43) 

so we have a double-counting of the available states of the system. We have four states of energy 0, and two of energy 
-E<j. To agree with physical results, we remove this double counting, so for each partition function written we divide 
the number of available states to the system by two. That is, for the single box case rather than Z = 4 + 2e^ Ed as 
one might expect from the degeneracies of the wavefunctions written, we will have Z = 2 + e^ Ed which agrees with 
our prior knowledge of this system. 

The rules for writing wavefunctions are fairly simple. Since 

[# (t) ,Q] = [0 (t) ,Y] =0^V|^> e{Q = Q ,Y = Y }36 ( V\v>e{Q = Q ,Y = Y } (44) 

[9^\X aa ,] = [i{00t _ e^e}, X aa ,} = 

^V|//> e{Q = Qo,Y = Y }BX aa ,\v>e{Q = Q ,Y = Y } (45) 

{» (t Uo,} = o [l{00t _ 6^},X 0a ] = 

=> V\v > e{Q = Q , Y = Y }3X 0r r \u >e {Q = Qo, Y = Y } (46) 

Where we have used the fact that the Hubbard operators conserve particle number and slave representation in the 
later two statements. This means that by finding a 'highest weight' state and applying these operators we can generate 
a basis of states consistent with Q = Qo, Y — Y 0} N = A . An additional subtlety to be aware of is that the Hubbard 
commutation relations mean that for blind application of this method we will end up with a non-orthogonal basis 
since: 

[ x n^ x o]] = -Xqi (47) 
so we will intentionally discard one of the three to form a basis. 



(c) Exact results at general N. To calculate the general partition function for the atomic model, we need to simply 
count the number of possible states of the system at each energy. This can be simply performed by writing the 
corresponding Young Tableau for each energy level and using the counting shown above in Eq. (8) (see Fig. 5). As 
one can only have one slave fermion (as it is a grassman variable), but up to h slave bosons, an L-shaped tableau 
of energy and dimensions (h,w) allows two L-shaped energy —Ed states to arise after application of the Hubbard 
operators with corresponding dimensions (h, w-1) and (h-1, w). Application of Eq. (8) yields the coefficient of the 
e P E d term. At energy — Ed we have Young Tableaux with dimension (h-1, w-1) and (h-2,w), and one can continue this 
procedure until energy —hEd when only one diagram (l,w-l) remains. This procedure allows us to generate exact 
results for the free energy in cases that would otherwise prove unaccessible to us. 
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FIG. 5. Illustration of counting procedure to generate the partition function of Q=3,Y=0. 



Using the above algorithm, it is simple to generate exact results for the free energy of an arbitrarily shaped Young 
Tableau (see Table II). 



Q,N 


F 


{3,4} 


-Tln(4) -Tln(5 + 4e' JEd +e 2 ^) 


{5,6} 


-Tln(21) - Tln(16 + 15e pE " + 6e 20Ed + e 3/3Ed ) 


{7,8} 


-Tln(120) - Tln(55 + 56e 0E " + 28e 2/3Ed + 8e 3/3Ed + e 40Ed ) 


{9,10} 


-Tln(715) - Tln(196 + 210e fiEd + \2Qe 2f3Ed + 45e 3/3Ed + We 40Ed + e 50Ed ) 


TABLE II. The Y=0 series with h=\. 


Using the above algorithm, it is simple to generate exact results for the free energy 


of an arbitrarily shaped Young Tableau. 
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2. Mean Field + Gaussian corrections 



(a) Susy. For the full atomic model, the (mean field) Hamiltonian reads 

H = -Edin^ + n x ) + X(n b + n f + n x + n<t> - Q ) + ((h f + n<p - (h b + h x ) + ^- < [6, <9 f ] > -Y ) (48) 

where n x and are slave fermion and slave boson fields respectively. As we saw in the free spin limit, it is useful to 
express Q and Y in terms of the height and width of the Young tableau (h,w), which slightly modifies the constraint 
on Q , as we absorb a term into the gaussian fluctuations. Then the free energy is 

F = F mft + 1f„. c (49) 

where 

F mft = -TIn(l + e -^ x "- E ^) - YTln(l + e~ px t) + Tln(l - e -0(A/-E-)) + NT ln(l - e~ 0Xb ) + T]n(e K + e~ K ) 
-X(h + w - 1) - £(h - to) 
= NT(h\n{n f ) + (1 - h) ln(l - h f ) +wln(n b ) - (1 + w) ln(l + n b )) 

+Tln(l-n x )-Tln(l + n^)-Tln( . n6n { ) (50) 



n b + n f 



and 



^F g . c = F 5Xb + F 5Xf 
T 

= +— ln((27r) 2 ((7Vn b (l + h b ) + n x (l - n x ) - n a (l - n a )){Nh f {\ -n f ) + n^(l + n ) - n a {\ - n a )) 
-nl{l-n a f)) (51) 

where n b = -g^- b — - and hf = ii^ s +1 ■ We have again dropped the Tin 2 from the double counting of the corner box 
and the familiar — Tln( ft "^/^ ) has again been included in the mean field to get us as close as possible to the correct 
saddle-point. The saddle-point solution d Q™* t = ^g^p = now yields: 

fif + jjUtp + ±(1 - n a ) = h (52) 
n b + jrn x + jrn a = to (53) 



where n a = e2f jl +1 = zr^i — gauges the degree by which a representation is fermionic or bosonic (refer to Appendix 



D for a discussion), = BEd( ± x) 1 and n x = -pE d , ■ Hence as in the free spin limit, we include as much 

as possible in the mean field saddle-point before treating the gaussian fluctuations. 

As in the free spin case, the saddle-point has an analytic solution for columns or rows. The special case Q =1, 
N=2 is of interest despite the non-physical divergence of the entropy noted in the free spin limit (Fig. 3). The heat 
capacity does not suffer from this malady, and is accessible to us as we can analytically (see Table III) remove the 
divergent terms in the entropy (which are independent of temperature-at least in the slave boson/slave fermion 
limits) by hand. Furthermore, the phase diagram generated by the saddle point solutions is interesting in its own 
right, with three possible solutions arising: slave bosons, slave fermions or mixed solutions. In Fig. 6 we sketch the 
simple picture arising. By minimizing the free energy (see Fig. 6(a)), one is driven to the cross-over from a slave 
fermion description to a slave boson description as one increases the value of Ed- This appears to be the physically 
realized situation. A more compelling phase diagram corresponds to maximizing the free energy of the saddle-point 
solutions, yielding a localized spin slave boson description, a small window of mixed solution followed by a slave 
fermion description at high temperatures which looks superficially similar to that seen close to quantum critical 
points (see Fig. 6(b)). Unfortunately it is not as simple to remove the non-physical divergence from this term, and it 
appears that simple procedures would lead to a non-physical heat capacity, such that the mixed solution likely is not 
chosen in the atomic model. It is interesting to note that at the mean field level the entropies match at these phase 
boundaries, and further that the difference between the entropy difference calculated from these approximations at 
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the level of gaussian fluctuations would lead to a weakly positive heat capacity, leading one to believe that if the 
gaussian fluctuations were properly controlled in the supersymmetric approximation, one might in fact see a second 
order transition. Unfortunately as yet we have not discovered a way to realize this goal, ft would be interesting 
however to imagine turning on a hybridization to a conduction sea and investigating if beyond some hybridization 
this state becomes realizable. For general Young Tableau of the row type we see that an analytic slave fermion solution 
exists, the column type has an analytic slave boson solution arising from the mean field consistency conditions. For 
more general L-shapes a numerical root-finding is necessary. 




FIG. 6. The analytic saddle-point solutions at N = 2 allow two possible phase diagrams: a) The physical saddle-point which 
picks the smallest free energy of the saddle-point solutions, b) A non-physical saddle-point which picks the highest free energy. 
While this looks qualitatively similar to that found close to a quantum critical point, we must remember that this model has 
no hybridization term. Non-physical divergences in a) are simple to remove analytically (constant in S) and do not affect 
the heat capacity, but this is less clear in the mixed phase of b)-simple attempts to remove this divergence seem to lead to a 
non-physically negative heat capacity. By turning on the hybridization, one might hope to garner a physical solution underlying 
this symmetry of the underlying formalism. 

















n f 


K . 


Hi 


n x 






n a 



i(l + e-f> E d) 2(l + e-P E d) U U 



2(l+e-' it: d) " " l+e-P^d 

2e-P E d-i l-e-f E d l-e-P E d 2 e -' 3E d-l 2(2 e -' 3E d -1) 

2{e-P E d+l) l + e-' 3E d l + e-f ,E d ( e -l 3E d + l) l + e -' 3E d 

TABLE III. Three solutions at N —2: slave fermion, slave boson, mixed. The latter only exists in the range Ed — to 
E d = Tln(2), has a constant ratio of slave partner to spin partner, and tunes continuously from slave fermion to slave boson. 
These analytic solutions allow us to obtain analytic solutions for the free energy, removing by hand non-physical divergences 
at N=2 for the slave boson and slave fermion limits. 
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(b) Slave fermions. For slave fermions, the Hamiltonian is 

H 8f = -E d h x + X((h x + n b ) - Q ) (54) 
such that the mean field free energy is: 

Fmft = -Tln(l + e-^ x ~ E ^) + NTqln{n b ) - NT(1 + q) ln(l + n b ) (55) 
and the gaussian corrections to this are 



^Fg.c = | ln(27r(iVn 6 (l + n b ) + n x (l - n x ))) (56) 



where n x = e -fiE d ^ + j^^ +1 is the number of slave fermions, and the particles have been constrained such that 

nb + -^n x = q (57) 

For this case, we see that as T — > for Ed < 0, n x — ► and n b — > g, while for £^ > 0, we have n x — ► 1 and 
fib ^ q — jj , n x (see Appendix D). 

(c) 57ave bosons. For slave bosons, the Hamiltonian is given by 

H sb = -Edhcj, + A(n/ +h ( j > - Qq) (58) 

which means that the mean field free energy is 

F mft = T ln(l - e -/»(*-®.0) + 7VT(1 - q) ln(l - n/) + NTqln(n f ) (59) 

and the gaussian corrections are 

jfFg.c = 1 111(27^/(1 - n f ) + n (l + n*))) (60) 
where = -0u d , _-L_n_i an( ^ * ne constraint maintained is: 

+ = 9 (61) 

For Ed > 0, as T — > 0, n/ — > and -^n^ — > while for £7</ < 0, — > and n/ — > Despite the complete 
transition to the bosonic vacuum state for E < 0, at low T some entropy does remain in the system. Both entropies 
remain positive in this approximation, while the entropy of the system with E > is much larger. 

B. Comparison of accuracy 

1. N=2, q = \ 

As these are all large N techniques, we note that this is an especially cruel limit. Nonetheless, a comparison shows 
that the heat capacity of slave bosons does the best job for Ed > (empty state favored) while the heat capacity 
for slave fermions does the best job for Ed < (spin state favored). Choosing the lowest free energy saddlcpoint 
corresponding to Fig 6 (a), the susy approximation agrees with the slave boson result for Ed > and slave fermion 
for Ed < at low temperatures (see Fig. 7). 
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C v vs. T for E d = 1 .5 (Empty) 




C v vs. T for E d = -1 .5 (Spin) 




The lowest free energy susy solution selects slave boson 




The lowest energy saddle point of susy is slave fermion 
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FIG. 7. (a) C„(T) for E d = 1.5, Q = 1, N=2;(b) C„(T) for E d = -1.5, Q = 1, N=2. When evaluated at the mean 
field saddle-point, susy is able to describe both the spin and charge dominated regions of the atomic model. A non-physical 
divergence can be removed by hand for the lowest energy saddle point solutions and does not affect the heat capacity. The 
potentially interesting mixed saddle-point solution does not have this nice feature and has been omitted as uncontrolled. 



2. The Y=0 series - results and errors 



In Fig. 8 we compare the heat capacity of our approximation with the exact result. We see that as we increase 
N while holding h = ^ and Y=0 constant the fit improves, although it is surprisingly good already at N=4. This is 
partly because the non-physical divergence of the susy approximation as N — > 2 does not affect the heat capacity as 
noted in section A2 (a). 



O vs. T for 



-1.5, 0=3, N=4, Y=0 



C„ vs. T for E . = 1 .5, 0=3, N=4, Y=0 



C u exact 

C u susy 




C„ vs. I tor 



-1.5, Q=19, N=20, Y=0 



C„ vs. T for E ri =1 .5, Q=1 9, N=20, Y=0 



C v exact 

C u susy 



C v exact 

C v susy 



FIG. 8. Heat capacity as N increases. 
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In Fig. 9 we see that the entropy being an extrinsic quantity improves dramatically in its description from N=3 to 
N=20 for the same Y=0 series. 



S vs. T for E d =-1 .5, Q=3, N=4, Y=0 



S exact 
S susy 



S vs. T for E d =1 .5, Q=3, N=4, Y=0 




S vs. T for E H =-1 .5, Q=1 9, N=20, Y=0 



S exact 
S susy 



S vs. T for E ri =1 .5, Q=1 9, N=20, Y=0 



S exact 

S susy 



T T 

FIG. 9. Entropy as a function of N. Since entropy scales with N, at large values it is approximately exact. 

But the eye is somewhat deceiving. In Fig. 10 and 11 we show the absolute errors the entropy and heat capacity 
respectively as a function of temperature. Both results seem to indicate that there is an error inherent to the 
approximation of the order of -t, which may also be present for the slave boson and slave fermion approximations 
and should be investigated. 
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FIG. 10. Absolute error in our entropy approximation as N is increased. (a) Ed < 0; (b) Ed > 
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Absolute Error in C v vs. T across N Absolute Error in C v vs. T across N 



E d = -1 .5, Y=0 E d = 1 .5, Y=0 




FIG. 11. Absolute error in our heat capacity approximation as N is increased. The sharp feature is likely a numerical problem 
where root-finding methods overlap. (a) Ed < 0; (b) Ed > 



4. CONCLUSIONS 

We have presented a simple strongly correlated electron model to test three representations of Hubbard operators. 
For the limiting case of spin fluctuations, we have demonstrated that mean field theory plus gaussian corrections is 
capable of capturing the leading order O(^) contributions for our model by comparison with Stirling's approximation 
expansions of exact results. Also in this case, by considering a mean field treatment plus gaussian corrections we 
have shown that a Schwinger boson description produces the most accurate results of the three large N methods. By 
introducing charge fluctuations, we have shown that this method of comparison can be extended to calculate heat 
capacities, and have found that spin states within this model generally are better described by bosonic descriptions 
(slave fermion) while empty favored by fermionic descriptions (slave boson), in the dangerous limit of N=2. The lowest 
energy saddle-point solution arising from the susy representation automatically makes this choice at the mean field 
level, but encounters a non-physical divergence in its gaussian fluctuations at N=2. The highest energy saddle-point 
of the susy representation on the other hand exhibits a qualitative similarity to that seen experimentally by systems 
near an antiferromagnetic quantum critical point, and it will be interesting to see if adding interactions to this simple 
model might drive us towards the realization of this type of phase diagram which might be expected to have two 
distinct bosonic order parameters: charge (j>Q which might couple to the hybridization and spin < b > which might 
allow the realization of a magnetically ordered state. Furthermore, we have shown that these results can be extended 
to higher N by the application of simple rules. We have seen that common to these approaches, the results become 
exact in the limit of large N, implying that if magnetism and paramagnetism co-exist in the ground-state of a material, 
at large N one might be justified in using the supersymmetric operators. 
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APPENDIX A: GAUSSIAN APPROXIMATION TO THE ENTROPY OF FREE SPINS 

The first term of Eqn (7) corresponds to the mean field theory in our approximate large N treatment. Since we 
cannot enforce the constraints Q — Qo and Y = Y"o rigorously at the mean field level, (Y has a term of order -^), 
we would like to see how inclusion of these terms consistently would affect the saddle point approximation we have 
made. The total free energy then is made up of five terms to order 

F tot = NFf + NF b + F n - XQ - (Y Q + F SXj + F SXb = F mft + (Al) 
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where the two terms of order N define the location of the large N saddle point, and the two latter terms are evaluated 
at the location of the mean field saddle point. The third term is the expectation value of the fermionic gauge field 9 
which we calculate at the level of gaussian fluctuations and then include in the saddle point, while fourth and fifth 
terms are corrections to the fermionic and bosonic free energies respectively. The first and second terms are simply 
the free fermionic and bosonic free energies: 

NFf = -NT ln(l + e~ l3Xf ) NF b = A/Tln(l - e _/3Ab ) (A2) 

The third term is the constraint term that we can not strictly handle at the mean field level, — ^ < [0, 9^] >= 
— + £• To accomplish this, we perform a Hubbard- Stratonovich transformation in the following way: 

^ + ^fj)(6 + 9Z V ) - |*t* = 9l m + + r>6 (A3) 

Since with the introduction of this grassman field we no longer have a quartic term in fields, we can treat these 
fluctuations at the gaussian level. This corresponds to calculating 




2( 2( 



U v = : B+ \_ r> =^7 + ^M (A4) 



where 

4>(iu; n ) = -T^G f (w + v)G h {y) = § ^ z+ ' (A5) 



_ n(-iu n +\ f ) n(\ b ) _ f{\f)+n(X b ) 



-iw„+2C iw„-2C UJ-2C w-2C 

Which means that 



(A6) 



_ Qo ,_Qo_ _ __Qou_ (A7] 

" 2C w-2C 2((uj-2() [ ' 
So the contribution to the free energy from these fluctuations is given by 

Fn = -J3 E ln(n(*w„)) = -i E ln(f ) + ± £ In(tw„ - 2C) - £ E M^n) (A8) 

= 0-f*&Hz-2C)+f*&M* + 0) (A9) 

= I ln(l + e- 2 ^ - I ln(2)) = I ln(i±^) (A10) 



Now, adding back in the additional £ from the constraint (as e Xn ^^) which was not done for the susy KondcEf 

model, this becomes F v = Tln(- — =y ). At the mean field level the constraint term reads — XQ = — A(n& + Uf — 1) 

which implies an extra contribution of A to the free energy hence this fluctuation term becomes F v = T \n( e>J f ^ e> b ) = 
Tln(^i). 

The corrections coming from varying the Lagrange multiplier about its mean field value and from the new gauge 

field arise from an expansion of the Lagrangian as Z = JV(f a ,b a ,5Xf 7 SXb,ri)e~-f dTCo (l + 5C + ^(SC) 2 ) where 
Lq represents the mean field Lagrangian, and we perform a saddle point evaluation. For a grassman gauge field, 
the expansion terminates at second order and the first order terms vanish so we have therefore e^ sc ^ , a grassman 
gaussian integral. For the Lagrange multiplier corrections, the first non-zero term is the second order diagram and in 
the gaussian approximation, we treat this as an exponential correction analogous to the fermionic case. This being 
the case, the corrections are as follows. The diagram of interest for Fg\ f is : 




(AH) 
This is given by: 
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T £ Gf{iuJn ) Gf{iUn) = Y J §^tjLf{z) J Eh (A12) 
= N{-ph f {\ s ){l-h f {\f))) (A13) 



where 



For a real gaussian variable, we have Z = f d ^ 5X} ^ e N fo c2t 0™/('V)( 1 n /( A /))( |5A /) _ 1 / j _ ~ — 

the integration is a sum over imaginary time (so it doesn't affect hf), which totals /3, the inverse temperature. The 
integration variable is normalized by (3 to make it dimensionless. The corresponding free energy contribution from 
these fluctuations is then: 

F SXj = -Tin 2 = |(ln(27riV7v(l - »/))) (AH) 

which is familiar as the gaussian fluctuations about the mean field for Abrikosov fermions. The diagram of interest 
for F S \. is: 




(A15) 



Which is given by: 



T £ G b {iv n )G b {i Vn ) = E f^n(^) I5 ^ F (A16) 

iv n ot—\..N a=l..N 

= E /^f^ i I^ = -A^ b (A fc )(l + n b (A b ))) (A17) 

a=l..JV 



Following the same procedure to find the partition function factor, we get Z = ^ y 2 Nn b (x b ){i+n b (\ b )) so * ne ^ ree 
energy contribution from these fluctuations is: 

F 5Xb = |(ln(27rAn b (l + n 6 ))) (A18) 

which is familiar as the gaussian fluctuations about the mean field for Schwinger bosons. Enforcing the constraint 
on Y at the mean field level, couples these two contributions as additional terms arise as 



dXjdXb dX b d\f dX* 



dX f 



F 

■f-=-prt a {l — n a ) leading to the combined contribution to the gaussian fluctuations as 



Nhf(l — hf) — n a (l — n a ) n a {\ — n a ) 

n a (l-n a ) Nh b (l + ii b ) - n a (l - n a ) 



(A19) 



The generalization of this term to the atomic model is straightforward as one incorporates fluctuations in the number 
of slave partners as well. 

APPENDIX B: Q=2, N=2 ATOMIC WAVEFUNCTIONS 

For Q = 2, N = 2, Y = 1, we have a totally antisymmetric wavefunction, so to write wavefunctions for these states 
we start with a simple two-particle anti-symmetric wavefunction and operate with lowering operators (Xo a , 0). The 
spin- spin Hubbard operator does not create new states because the initial wavefunctions are totally antisymmetric. 
We have with no slaves : 

|i >=////, e\A> (bi) 

whose wavefunctions associated are: 

/|/||0>, fe}//|0>-6{/ T t |0> (B2) 
With one vacuum state we have four possibilities: 

X 0L \A>, Y 0T |i>, eY oi |i>, 9Ao T |i> (B3) 
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/ T V|0>, /y|0>, (b\^ - x \fj)\0>, (6 t i 0t- x t / t ) | O> (B4) 
While with two empty states, 

A oi A OT |i>, ex ol x M> (B5) 

are the remaining possibilities: 

(0t)2| O >, 0V|O> (B6) 

For Q = 2, Y = -1, N = 2, we have a completely symmetric wavefunction, so we naturally start with a bosonic 
description as the highest weight state and apply boson lowering operators (Ao CT ,0^) and spin exchange operators 
(X aa ) to generate the complete set. The fully particle description is comprised of 

\B>=b\b\\0>, X n \B>, X n X n \B>,Q*\B>, Q*X n \B>, &X n X n \B> (B7) 

which in terms of fcrmionic and bosonic operators appears as: 

b\b\\0>, b\b\\0>, b\b\\0>, flb\\0>, / T t 6{|0>+/j6{|0>, $b\\0> (B8) 

Only fluctuations to one empty site are possible {Xq — 0)within the constraints for these quantum numbers. Thus 
the other possible states are categorized by: 

X ol \B>, A 0T |5>, X ol e^\B>, X ^X n X n \B> (B9) 
6 t i X t |0>, 6{x + |0 >, t 6 t i |O>+x t /||O>, t 6 t T |O>+ X t/t|o> (BIO) 

APPENDIX C: Q=3, N=2 ATOMIC WAVEFUNCTIONS 

For the (mixed) atomic model we can treat Q = 3 exactly (without comparison) for all three possible cases (Y = 
-2,0,2). Here for comparison's sake we list the wavefunctions of the N = 2 result, so that we can compare with a 
fairly simple calculation where the corner particle gets to choose its type in a non-trivial way. Since the degeneracy 
of the level gives us a value for the partition function, we can easily calculate a free energy and a heat capacity for 
this atomic model, and compare these results to the algorithm of section II i) d). For Q = 3,Y = 0,N = 2, our basis 
consists of four states with no slave particles, 8 = 2 3 distinct states with one slave particle and four states with two 
slave particles. For completeness I will list them here, and how they were derived. With no slave particles we have: 

|^>=&i/f/J|o>, e\A>, x u \a>, ex n \A> (ci) 

which are respectively: 

&{/ T VJ|0>, 6{(6}/J|0 > -&{/> >), 6 t T /|/||0>, 6{(6{/J|0 > -fcJ/flO >). (C2) 
With one slave particle: 

X 0] \A>, X ol \A>, OA 0T |A>, QX ol \A>, 

X 0] X n \A>, GX (n X n \A>, X ol X n \A>, QX ol X n \A> (C3) 

which are: 

<flb\fl X t /|/ i t |0> -^6{/f|0 >, 6l(^6{|0>-x + /J|0>), (C4) 

26{x + /f|0 > -6}(xVj|0 > +^6{|0 >), x + /f/J|0 > +0 + &}/J|O >, (C5) 
fel(x t / T t |O>+0 t 6 t T |O>)-26 t TX t/t|o>, -^6} /T i| >, 6}( x t/t|o>-0t 6 t| O>) . (C6) 

Here we have not written two equivalent eigenvectors. While with two slave particles: 

A j,A ||^4 >, QXoiX ]\A >, X^XoiX ]\A >, QX^XoiXo^ \A > (C7) 
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whose wavefunctions are: 

t (x t /||O>+^ t 6{|O>), 6|0t x t| O >, 0t( x t / t| O>+ 0t 5 t )j 6 t T 0t x t| O > (C8) 

We can also treat the fully symmetric case Q = 3, Y = -2, N = 2. As in the case with one less particle, we have 
spin fluctuations, wavefunction character fluctuations and limited charge fluctuations. The filled particle levels are 
described by: 

\C>= (6{) 3 |0 >, X tl \C>, X n X n \C>, X* l \0>, (C9) 

et|c>, Qix n \c>, e^x n x n \c>, etx T 3 jo> (cio) 

which are: 

(6{) 3 |0>, (6{) 2 fej|0>, (&{) 2 &{|0>,(&{) 3 |0> (Cll) 

/ i t (6{) 2 |0>, 6{(/f6{|0>+2/}6}|0>), 6} (/ t 6 t| 0>+2/ t 6 t|o>), /f(6}) 2 |0 > (C12) 

and with one less: 

X ol \C>, X ol X n \C>, X ol X^\C>, (C13) 

x ol ^\c>, x Ql e^x n \c>, x ol e^x^\c> (cu) 

X t (6 t i ) 2 |0>, x t b\b\\0>, xHb\f\0>, (C15) 
0t(&t)2| O>+2x t/t 6 t| O>) 0t 6 t 6 t| O>+x t (/ t 6 t| O>+/ t 6 || O>)) 0t( & t )2 | O > +2 x t/t 6 || > (C16) 

And of course we can treat the fully antisymmetric representation (Q=3,Y=2,N=2). Note that since N = 2 we can 
no longer have 3 spins! Here's what we have (totally unoccupied): 

\D >= (0t) 3 |O >, 6|D> (C17) 

which are 

(<^) 3 |0>, x t (0 t ) 2 |O> (C18) 

one spin 

X^ \D>, X l0 \D>, A T0 6| J D>, X lo e\D> (C19) 

which look like 

/ T t (0t)2| O >, /{(</»t) 2 |0 >, 6 t T (0t) 2 | O >+2/ T V ( / ) t|o>, 6t^t)2| 0>+2/ t x t^t|o> (C20) 

two spins: 

X^ X l0 \D>, X^X l0 O\D> (C21) 

and these are: 

/j/yiO >, <^(/» > -/» >) + xVj/flO > (C22) 
APPENDIX D: MEAN FIELD OCCUPATION 

To gain an appreciation of how the character of the spin system changes as a function of temperature, in Fig. 12 
we show the mean field results of the slave fermion approximation for Q = 1, N=2. We see (Fig. 12 (a)) that in this 
approximation when Ed > the Hubbard operators describe a charged hole which obeys grassman statistics. As the 
temperature is raised the bosonic description of the spin begins to carry more weight. When Ed < 0, the ground state 
at T=0 is a pure spin state which acquires a bosonic description (Schwinger bosons) (Fig. 12 (b)). 
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Slave Fermions: occupation vs. T (charge) 



Slave fermions: occupation vs. T (spin) 



Mean field values 



Mean field values 
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FIG. 12. Slave fermions: (a) E d = 1.5 (empty); (b) E d 



-1.5 



In Fig. 13 we show the mean field results of the slave boson approximation for Q=l, N=2. In Fig. 13 (a) we 
see that for Ed > the zero temperature ground state of the system is a positively charged bosonic hole. As the 
temperature is increased some fermionic spinorial character of the electron returns. When E<t < the system enters a 
fermionic ground state at low temperatures which is just the free spin case treated in section 2 (Abrikosov fermions). 
(see Fig. 13 (b)). 



Slave bosons: occupation vs. T (charge) 



Slave bosons: occupation vs. T (spin) 



Mean field values 



Mean field values 
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FIG. 13. Slave bosons: (a) E d = 1.5; (b) E d = -1.5 



For the case of physical interest (Q=1,N=2,Y=0), we saw that the supersymmetric description has three possible 
physical analytic solutions (see Table III). When n a — one recovers the slave fermion mean field results (Fig. 12), 
when n a = 1 one recovers the slave boson results (Fig. 13), and a mixed solution exists in the small region Ed = to 
Ed = Tln(2) as shown in the inset to Fig. 14. If one follows the Y=0 series by increasing N while holding h = = i 
constant (see Fig. 14), one sees that one generically has a mixed solution - neither slave bosons nor slave fermions can 
describe the physics of systems exhibiting this symmetry. As approximated methods become asymptotically exact as 
N — > oo, if the phase diagram of Fig. 6 (b) were to become the physically realized state, one might be able to justify 
a large N expansion to treat the physics of the mixed region using susy Hubbard operators at Y=0. 



22 



0.8 



0.6 



0.4 



0.2 



Y = Series: n a vs. T 

Can we access a region hinted at in N=2? (N, E d ] 
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FIG. 14. The Y = series. In the free spin limit (Ed/T « 0), the corner box favors a fermionic description of the spin 
(slave bosons n a = 1) while in the free charge limit (Ed/T >> 0), the corner box favors a fermionic description of the charge 
(slave fermions n a = ). At high temperatures a mixed solution exists slightly favoring slave bosons. Both the heat capacity 
and entropy approximations for N>2 produce physical results and may be seen as an extension of the high temperature N=2 
mixed solution (see inset) which for Ed > and T > j^|r tunes from slave fermion to slave boson as temperature increases. 
Study of the physics of this series in less trivial models may reproduce some of the interesting physics seen close to a quantum 
critical point. 
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